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Abstract 

The Chiral Random Matrix Model or the Gaussian Penner Model (gen- 
eralized Laguerre ensemble) is re-examined in the light of the results which 
have been found in double well matrix models [[l], |2[] and subtleties discovered 
in the single well matrix models ||. The orthogonal polynomial method is 
used to extend the universality to include non-polynomial potentials. The 
new asymptotic ansatz is derived (different from Szego's result) using saddle 
point techniques. The density-density correlators are the same as that found 
for the double well models ref. |2j (there the results have been derived for 
arbitrary potentials). In the smoothed large N limit they are sensitive to 
odd and even N where N is the size of the matrix |2[]. This is a more realistic 
random matrix model of mesoscopic systems with density of eigenvalues with 
gaps. The eigenvalues see a brick-wall potential at the origin. This would 
correspond to sharp edges in a real mesoscopic system or a reflecting bound- 
ary. Hence the results for the two-point density- density correlation function 
may be useful in finding one eigenvalue effects in experiments in mesoscopic 



systems or small metallic grains. These results may also be relevant for 
studies of structural glasses as described in ref. ||. 
PACS: 02.70.Ns, 61.20.Lc, 61.43.Fs 
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1 Introduction 



In this work we re-visit the Penner model in the context of disordered meso- 
scopic systems where they naturally appear e.g. in disordered models of 
metal-insulator transitions |6|, [7|, superconducting-normal interface with An- 
dreev reflection (the openning of a gap) |§ and also as models of structural 
glasses 0. Historically the Penner matrix model was studied in the context 
of the moduli space of a punctured surface 0. There an equality between 
the Penner matrix model and the Euler characteristic of moduli space of 
punctured surfaces was computed, before taking the continuum limit. Later 
in [irOJ it was shown that even after taking the continuum limit the Euler 
characteristic of moduli space of unpunctured surfaces was obtain as the 
free energy of the Penner model. This was done in an effort to understand 
2-dimensional gravity coupled to matter at the critical point c=l. Recent ad- 



vances in this direction are contained in [111. There has also been recent work 



on the generalized Laguerre ensemble in the context of the Chiral Random 
Matrix models of QCD (see [12, 13] and references therein) and in describing 
a novel group structure associated with scattering in disordered mesoscopic 
wires @. 

This study focuses on the Gaussian Penner model with potential V(M) = 
\liM 2 - f InM 2 where M is a N x N random matrix; after some work on 
gapped matrix models has been reported and clarified ||. In the recent 
work two basic ideas have been implemented. First the idea of symmetry 
breaking (see || £§]) and then a new asymptotic ansatz of the orthogonal 
polynomials Here for the Gaussian Penner model the corresponding poly- 
nomials are the associated Laguerre polynomials, L%(x), where both the a 
and N asymptotes are to be taken, in previous work only the N asymptote 
was found. It turns out that it is in this asymtotic region that the singular 
models make contact with the double well random matrix models and the uni- 
versality of the gapped matrix model is extended to include non-polynomial 
potentials. Ideas of symmetry breaking are also true for the Gaussian Pen- 
ner model 0, U but because of the brick wall potential at the origin, which 
involves one electron quantum tunneling through an infinite barrier, it is 
harder to give an intuitive picture. These ideas will be ellaborated on in 
future papers. 

The paper is divided as follows. It starts by describing completely the 
model and establishing the notation and conventions. In section 2, the old 
asymptotic ansatz of Szego for the associative Laguerre polynomials is taken 
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Figure 1: The potential for the Gaussian Penner Random Matrix Model with 
V(M) = \^M 2 - |lnM 2 , // = 1 and t = 1. 



and shown to not correspond to the new universality of the double well matrix 
model studied in |], @. In section 3, the new asymptotic ansatz is derived 
using saddle point techniques. This corresponds to the asymptotic ansatz 
discussed in [0, @. Section 4 ends with conclusions and open questions. 



2 The Model, Notation and Conventions 



We consider models of the type (see [14] for details of notation and defini- 
tions) 

Z = exp(-F) = J dMe- NTrV M ( 2 .1) 

where V{M) = V (M) - § InM 2 and here V (M) = ±/iM 2 . See Fig. |l[ 
In general the orthogonal polynomials are P n (A) = A n + l.o. where A are 
the eigenvalues of M and the orthogonality conditions are fj^ G?Ae~ w ( A )P n (A)P m (A) 
h n 5 nm . The partition function can be expressed in terms of the h n 's as 

Z = N\hohih 2 h N ^i. 

For the large N limit, the density of eigenvalues, p(z) = (-^) X)i=i 8(z— Aj), 
can be found by solving either a saddle point equation or the Schwinger- 
Dyson equation. In terms of the generating function F(z) = ■h 
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Figure 2: The density of eigenvalues for the Gaussian Penner Random Matrix 
Model for \i = 1 and t = 1. 



/ dz'^^f the Schwinger-Dyson equation reads F(z) 2 — V'{z)F{z) = M(z), 
where M(z) is a meromorphic function. The density of eigenvalues is given 
by p(z) = —-lmF(z). The generating function F(z) in the large N limit for 
t > is 

(2.2) 



where a 2 = M + ^(1 + 1) and 6 2 = 
corresponding density of eigenvalues. 



(2+£) 



(1 + t). See Fig. |2| for the 



3 Exact Solution for the Symmetric Gaussian 
Model 

The orthogonal polynomials satisfy the recurrence relation 

(Z - S n )P n ( Z ) = P n+1 ( Z ) + RnPn-^z) (3.1) 

where S n , R n are the recurrence coefficients. For symmetric models S n = 0. 
For even potentials instead of recurrence relation 

zP n ( Z ) = P n+1 ( Z ) + RnPn-^z) (3.2) 
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one can use 



Z 2 P 2n (z) = P2n+2{Z) + (R 2n +1 + R 2 n)P 2 n{z) + R 2n - X R 2n P 2n - 2 (z) (3.3) 

and 

Z 2 P 2n +l(z) = P2n+3(z) + (R2n+l+R2n+2)P2n+l(z)+R2n+lR2nP2n-l(z)- (3.4) 

Where we have multiplied z to Eq. (|3.2|) and expanded. Then eq. ( |3.3| ) 
contains only even polynomials P 2n (— z) = P 2n (z) and eq. ( |3.4| ) contains only 
odd polynomials P 2n+ i(—z) = — P 2n +i i z ) ■ This simplifies the solution as we 
will see. 

(1) . Let us first work with the even set. Let y = z 2 and define functions 
Pn(y) = P2n(z). In terms of these 'new' polynomials 

yVn(y) = V n +i(y) + S n V n {y) + TZ n V n -i(y) (3.5) 

where S n = R 2n + R 2n +i = A n + B n and lZ n = R 2 n-iR 2 n = A n B n _i. 
These polynomials obey 

/ dye- N '^~ t,l °^V n {y)V m {y) = h n 5 n , m (3.6) 
Jo 

where V (y) = 2V (z) = fiy + , t' = t - ^ and N' = f . 

This is the same as the brick wall problem i.e. the linear Penner model 
if t <-> If, N <-> N' and t'N' = (t - ^)f = ^f^. 

(2) . A similar analysis can be carried out for the odd set. Define 

V n (y) = z- l P 2n+1 (z). (3.7) 

Then 

yPn(y) = Pn+i(y) + s n v n (y) + n n P n ^{y) (3.8) 

where 

S n = R2n+l + R 2 n+2 (3-9) 

and 

T^n — R2n+\R 2 n- (3.10) 
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Because of the extra factor of z associated with the odd series the 'barred' 
polynomials satisfy orthogonality condition 

" dye~ N '^- Fl ^V n (y)V m (y) = K5 n , m (3.11) 

where t' = t + ^r- This barred system can be solved as a brick- wall 
problem. Original recurrence coefficients are obtained by 



R 2n+1 = -{S n + ^JSl - ATZ n } (3.12) 



R %n = -{S n -sJSl-4K n }. (3.13) 

Using W n+1 +W n +S n y n = ia±j±m an d S n W n+1 - W n - £ = -R n+1 y n+1 + 
R n yn-i (see ref. |TJJ for details and definitions of W and y) we get 

2n + l + t'N' An + l + tN 

s - = w = (3 ' 14) 

//iV' //TV v ; 

_ nfo + t'iV') _ 2n(2n-l + tN) 

n(n + W) _ 2n(2n + l + tAQ 

(1). For the even set the orthogonality relation 

dye- N '^- t/lo ^V n (y)V m (y) = K5 n , m (3.18) 

JO 

simplifies to 

/ dye- NVo <»)y»'«p n (y)p m (y) = h n 5 n , m . (3.19) 
Jo 

Now as 

n(n + t'iV') 2n(2n-l + tAQ 

^2^2 ~ ^2 N 2 > ^- 2U ) 
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_ 2(2-l + tAQ 2(1 + tN) 

etc. Note the following 

U n = (3.23) 



Thus h n is 



Note: 



for n = 



Then 



For 



/i n — 1Z n h n -\ — 1Z n 1Z n -ih n -2 

n(n + UN') (n - l)(n - 1 + t'iV') 
^ 2 7V' 2 ^iV' 2 

(1 + t'TV 7 ) 



/i 

4-i'M' -i- 1\ 

ho. (3.24) 



^ 2 iV' 2 

ra! r(n + t'iV + l) 



(/x7V') 2n T(t'N' + 1) 



Pn(j/) = P2n(z) (3.25) 

7> n (y) = Po(^) = 1. (3.26) 

/'OC 

= / dye- N ' v ^y N,t 'Vo(y)Vo(y) 
Jo 

dye- N ' v °Wy N,t '. (3.27) 

V (y) = HV, (3.28) 
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Defining 



/'DC 

ho = / dye- N '™y N,t ' . (3.29) 

JO 

u = N'fiy, (3.30) 

rfw = N'fxdy (3.31) 



we get 



k ° ~ Jo N^ 6 W 



7^FTi r (N't> + 1) (3.32) 



and 



(N'/jl) 



n\T(n + l + t'N>) 

(3.34) 



Rewritting 



V n {y) = Q n {u) (3.35) 
the orthogonality condition is 



/ duu N,t 'e- u Q n (u)Q m (u) = (N'fx) N,t ' +1 hJ n , m 
Jo 

(n\) 2 T(n+l+t'N') 



On,r. 



(fiN') 2n T(n+l) 

Redefine Q' n (u) = (j ^Q n (u). Then 

T( n+ 1 + t'N ') 
T(n+ 1) 

comparing with the Associated Laguerre Polynomials L^(x) we get 



f duu N,t ' e- u Q' n {u)Q' n {u) = ^ + 1 + (3.36) 
Jo 1 (n + 1 
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then (a = N't' Associated Laguerre Polynomials) 



V n {y) = Qn(u) = 



(fiN'y 



and 



(3.37) 



(3.38) 



My) = 



\/K. 



(fiN') 



AiV't'+l 



-r(n!)(-irLf(JV M y). (3.39) 



\| (n!)r(n + 1 + W) 
Thus the normalized even set of orthogonal polynomials are 



_i — J gy~5~e~g~L„ * (iV uy). (3.40) 



(2). For odd set 



i _ _ 

dye- N '^- t,l °^P n (y)P m (y) = h n 8 n 



(3.41) 



where t' — t + ^ and N't' = ^f^- The orthogonality condition are 



roo _ _ _ 

/ dye~ N ™y N 1 P n {y)P m {y) = h n S nm . (3.42) 
Jo 



Note that 



h 



For 



we get 



h n -i 
n(n + i') 



(3.43) 



(3.44) 
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hn T^nhn—1 

= 1Z n 1Z n -i1Z n -2 T^iho 

n\(n + i'N')\ - 
(fiN') 2n {fN')\ °" 

Note: 

P n (y) = z^Pun+i^) and P (y) = z~ x P x {z). 

But Pi (z) = z hence Pq(?/) = z ~ lz — 1- Therefore 



Defining 



and 



we get 



/in 



dye- N '^y N ' p Vo(y)V (y) 



poo 
JO 



/i 



It = iV'/ig/ 



du = N'fidy, 



r du -«/ " \N'F 

Jo N'n y N'fi } 

1 ^^e-^'*'. 



(N' fi) N ' t,+1 Jo 



So that 



/ir, 



r(n + l)r(n + 1 + W) 



(A^» 2 "+ JV ' f '+ 1 
Now let Qn('u) = V n (y). Therefore 
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poo _ _ 

/ due- u u N,t 'Q n (u)Q m (u) = (N'fi) N ' t,+1 h n 5 nm 
Jo 



r(n + l)r(n+l+t'N') . 

Or, 



(jiN') 



2n 



n! T(n+l+i'N') 
= (p^T) ~ { ^-(3.51) 



Define Q' n (u) = { -^Q n (u) then 



f° due- u u NT Q' n (u)Q' m (u) = {U+ , + ^ nm - (3.52) 

Jo n! 



For 



QM = 



where a = N't'. Now 



hence 



Vn{V) = Qn{u) 



TV. 



-Qn{u) 



:-i) n if(« 



n\{-l) n TN ' V/ATl , 



(jiN') 



(3.53) 



(3.54) 



P n (y) 



P n (y) 



\r(n + l + t'N>) { l) Ln {Nfiy) - 
Thus the normalized odd orthogonal polynomials are 



(3.55) 



= [ 



^..^ L s, {N , fl!jy (3 . 56) 



T(n + 1 + N't') 



—-\2y 2 e 2 
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The important lesson to learn from this exercise is that a = < - jVt ~^~ 1 - > — - 
for n' an integer which is even or odd and x = N'ny for the orthogonal 
polynomials, which are proportional to the generalized Laguerre polynomials 
(x) , of the Gaussian Penner Matrix Model. 

4 The Old Asymtotic Formula for Associa- 
tive Laguerre Ensemble 

Consider the asymptotic formula given in Szego's book on orthogonal poly- 
nomials. The formulas of Plancherel-Rotach type for Laguerre polynomials 
for a arbitrary and real, e fixed positive number for x = (An + 2a + 2) cos 2 0, 
e < < f — en~ are given below 



e 2 X L^\x) = (—l) n (nsin(f)) ? x % * 
{Cos[{n + ^-^){sin24> - 20) + j] + (nx)^0(l)}. (4.1) 

For the Gaussian Penner model with a = ( - Jvt ~^~ 1 ' > — - and x = N'fj,y 



e ^ L^\x) = (— 1)" (ir sirup) * x 2 ° ^n" a 

I ( Af-f _ (-l) n ') 7T 

{Cos[(n + -)(sm20 - 20) + ^ 1 — ^(sm20 - 20) + -] 

+(nx)^0(l)}. (4.2) 

The combined asymptotic ansatz for the even and odd results eq. ( |3.40| ) 
and eq. ( |3.56|) for the Gaussian Penner model in this asymtotic region is (n' 



stands for both even 2n and odd 2n + 1 integers) 

i>n'{ x ) — e ~^L^\x) = (— l) n (7rsm0)^ i x^~^n^~2 

{Cos[(n + - + ^)(sm20 - 20) + | - t^— (sm20 - 20)] 

+ (nxy^O(l)} (4.3) 

Note that the term ^— j— (sm20 — 20) is not the extra term (with rj see 
ref. |L], |2|] or eq. ( |5.15| ) below) that was found for the asymptotic ansatz in 
the double well problem. 
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5 The New Asymptotic Formula for Associa- 
tive Laguerre Ensemble 

Let us start with an expression which has a well defined a = limit. Using 

jL /(l) = i. /(u + x)U (5 ,) 
we can write using Cauchy's theorem 



X J 

nl au 

= ^(^) n [(^ + «) {n+Q) e- (u) ]|„=o 
n\ au 

Jc 2t7T Z n+1 ^ ' 



(5.2) 

in which the countour C is a small circle around the origin {\z\ < x) 
change z — > nz (\z\ < ^) 



f(z) = z + logz - (1 + -)log(z + -) (5.4) 

n n 

we explore - and - finite. In this representation the limit a = is well 

1 n n 1 

defined. C is a circle around the origin. The saddle point is given by an 
expression 



z 2 + z( ) + - = 0. (5.5) 

n n 



With the parametrization 



= 2W-cos0 (5.6) 

n V n 
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z = f-e* (5.7) 

and z~o are the saddle points. This is valid in the range \ SL ^\ < 2y^; 
if the parameters are such that this inequality is reversed, there is only one 
real saddle-point. This is where the new saddle points will have to be taken 
into account. 

If there is a saddle-point in the integral 

^-I e -»/W (5.8) 
2^7r z y ' 

at z then 

, e -n/(^0) / 2^ 

e 2 (5.9) 



zo \n\f"(z )\ 
with f"(zo) = \f"(zo)\e l6 . Indeed expand f(z) around z 

f(z) = f(z )+ 1 -(z-z o r\f"(z )\e^ (5.10) 
path is z — z = xe 2 

/ « e -*f-»/(«o) [ + °° dxe~^ f "^. (5.11) 
j —00 

Here we need to add the contributions of z , z and the other saddle 
points. Therefore if f(z ) = A + iB 



with 



e_ "V^b 2cos( " B+ ^ + ^ )(i) (5 ' i2) 



, fn2i sirup e ** 

f (zo) = \ 7= — (5.13) 

V x /£. 4- e}4> 



and 9 is it's phase. The (—1)™ would only show up in B. This is a new 
asymptotic expansion of the generalized Laguerre ensemble in a novel asymp- 
totic regime. It would be nice to have a physical picture of this asymptotic 
regime. 
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On simplifying the above expression it will give the asymtotic ansatz 
found in ref.[|l] for the double well matrix model. Following ref. Jl| for iV 
large but N — n « 0(1) and x lying in the two cuts the asymptotic ansatz 
for the orthogonal polynomials of the Gaussain Penner Matrix Model can be 
approximated by 



lfj n (x) = -jj 



cos(iVC - (N - n)<j> + X + (-lTv)(x) + 0(±) 



(5.14) 



where /, £, <f), \ and r\ are functions of x and are given by 



s n (b 2 — a 2 ) , . . 

f(x) = — 1 2 ; sin20(x) 



C'( x ) = -irp(x) 
cos20(x) = 



2 _ (a 2 +b 2 ) 
2 

(fe 2 -q 2 ) 



2 




COS 01 


X) 


X 




sin 


x) 



cos 277(0;) 

sin2^(x) = a 

x 

1 7T 

with a 2 and b 2 as given in sec. (2) for the Gaussian Penner model. 

All the corresponding correlation functions of this model will be as obained 
in ref. QXJ] and ref. . Following ref. and using the contour of integration 
for the Gaussian Penner model, see Fig. |3], the the smoothed density-density 
correlation function can be derived in the thermodynamic limit and is an 
oscillating function of A^: 

2vr 2 A^(A,^) r ^< 1 



Xfi(\fi -a 2 - b 2 ) + a 2 b 2 + (-l) N ab(fi - A) 2 ) . 

(5.16) 

Here for the symmetric Gaussian Penner Model a(z) = (z 2 — a 2 ){z 2 — b 2 ), 
a 2 = ^ + I v^TTT), b 2 = fi±a - l ^TT), ex = +1 for b < X < a, 
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Figure 3: The contour of integration for the two-cut Gaussian Penner model. 
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e\ = — 1 for —a < X < — b and j3 = 1,2,4 depending on whether M the 
matrix is real orthogonal, hermitian or self-dual quartonian. 



Using this expression in the formulas for the mesoscopic fluctuation flo 
would give rise to terms depending on N. This may be observed in single 
electron experiments on mesoscopic samples which have gaps in their eigen- 
value spectrum. Work to explicitly obtain all these expressions in the Chiral 
Matrix Model in this asymptotic regime is in progress. 



6 Conclusion 

Finally let us note that in matrix models when the number of connected 
components for the support of the eigenvalues changes, one finds a new uni- 
versality class for the correlators which has been extended here to include the 
non-polynomial potentials. It is not completely obvious that it is legitimate 
to use the simple one cut function in the application to mesoscopic fluctua- 
tions as the correlation function are different for these gapped random matrix 
models because of single electron tunneling effects. Further questions of sym- 
metry breaking in these singular random matrix models are open questions 
to be explored in the future. 
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